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Surface plasmon polaritons scattering from 2D subwavelength silicon wires is investigated using 
finite difference time domain method. It is shown that coupling an incident surface plasmon polariton 
to inter-cavity modes of the particle can dramatically changes transmitted fields and plasmon- 
induced forces. In particular, both transmission and optical forces are highly sensitive to the particle 
size that is related to the excitation of whispering gallery modes or standing-wave modes depending 
on the particle shape and size. This features might have potential sensing applications. 


I. INTRODUCTION 

Surface plasmon polaritons (SPP’s) are evanescently 
confined modes that propagate at the interface between 
a dielectric and a conductor. This confinement leads to 
an enhancement of electromagnetic fields at the interface 
which is very sensitive to the surface defects [lj |2]. The 
unique properties of the SPPs give them two dimensional 
nature and opens an opportunity for scaling down opti¬ 
cal devices to nanometric dimensions In order to 

understand basic behaviors of such devices the interac¬ 
tion of SPPs with small particles and defects is a problem 
of fundamental importance. While, scattering by surface 
defects in all-metal structures has been studied in many 
works mm, scattering by dielectric objects placed 
near a metallic surface has been explored less so far. The 
study of the latter structure might be particularly useful 
in the plasmonic optical tweezers in which a small object 
is trapped near a metallic surface [MO]' spp scatter¬ 
ing by finite-size gold nanocubes placed in the vicinity 
of a metal surface has been analyzed in Ref. ED. and 
scattering by 2D rectangular dielectric particles has also 
been analyzed |5]. 

While silicon is widely used in conventional photonic 
devices, it also provides an important material for the 
design of silicon plasmonic devices with advanced func¬ 
tionalities. Resonant light interactions of high-refractive 
index dielectric nanoparticles have already been analyzed 
in details nans]. Here, we consider their interactions 
with surface plasmon polaritons. It is well-known from 
Mie theory that for dielectric particles with higher refrac¬ 
tive indexes, the light scattering efficiency increases. Due 
to these properties, high index dielectric material offer 
great practical advantages in many plasmonic structures 
such as low loss propagation and strong field confinement 
in the hybrid plasmonic waveguides nuna and enhanced 
light transmission and enhanced trapping forces in the 
plasmonic nanoaperture traps mm- 

Dielectric microcavities can support whispering gallery 
modes (WGMs) which have been intensively applied to 
many devices such as optical sensors m , narrowband 
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filters Esa. and microlasers ED- Most research on WGM 
cavities are restricted to large cavities with sizes much 
larger than light wavelength. However, to develop ad¬ 
vance nanophotonic elements, it is of great importance 
to consider nanosize resonators with sizes smaller than 
the considered light wavelength. In this regime there are 
few research on WGMs (see, for instance [22/23;). In this 
paper we consider subwavelength silicon wires and inves¬ 
tigate the excitation of WGMs with small mode numbers. 

From the theoretical perspective, the scattering of sur¬ 
face plasmon polaritons by surface defects has been stud¬ 
ied by several techniques such as the Green’s Function 
method mm and Rayleigh expansion method m EiEa. 
On the other hand, the finite difference time domain 
(FDTD) method is a powerful and well known technique 
to determine electromagnetic fields around nanostruc¬ 
tures of complex and arbitrary geometries p6h29j . In this 
paper, we use two-dimensional FDTD method (a home¬ 
made code) for analyzing SPPs scattering from subwave¬ 
length silicon wires with different cross sections. We also 
use the Maxwell stress tensor (MST) technique to inves¬ 
tigate the optical forces on the particles. We show that 
SPP scattering parameters and optical forces on a sub¬ 
wavelength wire are non-trivial and strongly depend on 
its cross section shape and size. Indeed, our calculations 
show that the excitation of intercavity modes such as 
WGMs and standing wave modes, is responsible for this 
behavior. It is important to note that the high refractive 
index plays a key role here and for sufficiently high re¬ 
fractive index material such as silicon, the subwavelength 
wire would exhibit WGMs (with small mode numbers) 
even if its size is smaller than the vacuum wavelength of 
the light. Although they are less investigated, such low 
order modes in subwavelength regime could be of great 
importance for nanophotonics and nanosized devices [22]. 


II. CALCULATION OF OPTICAL FORCES 

To excite SPPs in the FDTD simulations we use the 
total-field/scatter-field technique [26]. A small grid size 
(less than 3nm) is chosen in order to account for all 
plasmonic near-field behaviors. We bound the simula¬ 
tion domain with a convolutional perfectly matched layer 
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(CPML) [26]. For the dielectric constant of the metal¬ 
lic region, we use the Drude model with parameters that 
fitted to the experimental data of the complex dielectric 
constant for silver [30] . 

Once the electromagnetic fields have been calculated 
using FDTD method, we can calculate the optical forces 
acting on the particle using the Maxwell stress tensor 
(MST) method [31] . The time-averaged force acting on 
the center of mass of the particle is 

(F) = f (T) • n dS (1) 

Js 

where S is a surface enclosing the particle, n is the unit 
vector perpendicular to the surface, and (T) is the time- 
averaged Maxwell stress tensor for harmonic fields, i.e., 

(T) = ^ jeE E* + M H H* — * (e|E | 2 + M |H| 2 ) J , 

( 2 ) 

where E is the electric field, H is the magnetic field, 
and e and p are the permittivity and permeability of the 
surrounding medium. 

We notice that for small particles or small entities in 
the structure (like small gaps) and for low refractive in¬ 
dex discrepancy between dielectric particle and surround¬ 
ing, the force calculations using the Maxwell stress ten¬ 
sor need a lot of care. In order to ensure the correct¬ 
ness of our results we verified the numerical convergence 
of our algorithms. In our simulations we considered a 
uniform two-dimensional square-cell space lattice with 
Ax = Az = A. For each configuration we performed 
a convergence analysis to ensure that the mesh size is 
chosen properly (i.e. it is small enough) to reduce nu¬ 
merical errors. As a result of this analysis we decided to 
employ a mesh size between 2 nm and 3 nm. 


III. COMPARISON OF FDTD COMPUTATIONS 
AND RAYLEIGH EXPANSION 
APPROXIMATION 

It is instructive to first compare FDTD computations 
with a theoretical model. We consider the scattering of 
SPPs by small defects on a metal-vacuum interface. As 
shown in Fig. la SPP of frequency co propagates along 
the x axis. Suppose that the metal-vacuum interface is 
describe by the following function: 


z = h(x). (3) 

This form represents a one-dimensional defect in a single 
interface all-metal structure. For small defects, the mag¬ 
netic field in the vacuum half-space can be written as a 
sum of the incident field and scattered field according to 
the conventional Rayleigh expansion [32 ) : 

H y (x , z) = exp[ik p x — k^z\ + J dkH y (k)exp[ikx — k z z\ 

( 4 ) 
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FIG. 1. (a) Schematic of the surface plasmon scattering from 
surface defects. Transmittance as a function of defect width 
for Gaussian (b) and rectangular (c) shape indentations in 
the silver-air surface (A — 600 nm and hmax — 3 6nm). Solid 
curves correspond to Rayleigh approximation and circles rep¬ 
resent the corresponding transmittance computed by FDTD 
method. 


where k p = y yqy, i+i are th e wave number 

of a surface plasmon polariton at a flat metal-vacuum 
interface and the inverse decay length of the fields per¬ 
pendicular to the interface, respectively. Also, e is the 
dielectric function of the metallic half space, ko = u)/c, 


and k z is defined by k z = ^k 2 — fcg, Re{k z } > 0. 

In order to find scattered coefficients, H y (k ), assuming 
that |e| 1, one can apply the surface impedance bound¬ 

ary condition (SIBC) to the tangential components of the 
fields at metal surface 0. As a result, the reflected and 
transmitted SPP fields may be written in the following 
form 


H y {x —»• oo, 0) = (1 + r)exp(ik p x) (5) 

Hy{ x —>> — oo, 0) = exp(ik p x ) + pexp(—ik p x ) (6) 

where r and p are parameters to be calculated using the 
boundary conditions [6 . Therefore, the transmission and 
reflection coefficients are 

T=|l + r| 2 , R=\p \ 2 (7) 

Also, the energy conservation law leads to scattering co¬ 
efficient from S = 1 — T — R. 

Figure [l] illustrates the results for two defect of rectan¬ 
gular (b) and Gaussian (c) profiles in a silver surface (at a 
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FIG. 2. (a) Transmittance and (b) reflectance of a circular 
silicon cylinder as a function of its diameter near a flat silver 
surface at A = 600 nm. (c)-(g) near magnetic field patterns, 
\H y |, in z > 0 half-space for the particles are shown by arrows 
in Fig. 2 (a) . The arrows show the Poynting vectors. The 
gap between the particle and silver surface is 12 nm. 


wavelength of 600 nm). We also perform two-dimensional 
FDTD method to compare the results. As seen in the fig¬ 
ure, there is a noticeable discrepancy between the FDTD 
results and the theoretical prediction of Rayleigh approx¬ 
imation. This discrepancy might be attributed to the 
simplified boundary conditions and ignoring the propa¬ 
gation losses in (§. Furthermore, here the defect depths , 
hmaxi is relatively large. Note that, for rectangular defect 
there is much less agreement which could be attributed to 
the violation of the surface impedance boundary condi¬ 
tion at the discontinuity points. However, for smooth and 
sufficiently shallow defects one would expect the Rayleigh 
approximation leads to more accurate results [6]. 


IV. SCATTERING BY SUBWAVELENGTH 
PARTICLES 


Let us now apply FDTD method to study the scat¬ 
tering of surface plasmon polaritons by subwavelength 
silicon wires. We first consider scattering of a SPP from 
a circular silicon cylinder placed near a flat silver surface. 
The surrounding medium in the z > 0 half-space is water 
(refractive index n m = 1.33), and the refractive index of 
Si was taken from Ref. [33] . We fixed the gap between 
the particles and metal surface as 12 nm. Figs. [5] (a) and 
(b) show transmittance and reflectance as a function of 
the particle diameter. As can be seen the particle behave 
in two different ways depending to its diameter. First, 
for very small particles (D /A < .15) most incident SPP is 


transmitted, however larger particles transmit less power. 
Second, for larger particles ( D/X > .15) the transmit¬ 
tance and reflectance curves strongly fluctuates. In this 
case, the transmittance suddenly rises and falls and sur¬ 
prisingly, in some points it is very high and in some points 
it is approximately zero. To better understand these be¬ 
haviors, the magnetic field distribution (| H y |) and Poynt¬ 
ing vector field are shown in Figs .[21(c) -(g) for some par¬ 
ticle of interest (depicted in Fig. [2](a) with arrows). For 
instance, Fig. [2] (c) shows excitation of second-order (az¬ 
imuthal) transverse magnetic (TM) polarized Whisper¬ 
ing Gallery mode (WGM) by evanescent coupling of in¬ 
cident SPP (here, the radius of the wire is 130 nm, which 
corresponds to Diameter/wavelength=0.433). More pre¬ 
cisely, the whispering gallery modes, ? </>), of a 2D 

cylindrical cavity of radius a can be analytically found by 
solving Helmholtz equation using the separation variable 
technique [34]. The solutions take the form, 

ipN,m{r, 4>) = A N>m e ±zm,p G m (k N:m r) , ( 8 ) 

where the function G m is Bessel function of first kind, 
J m , for r < a and the Hankel function of the first kind, 
Hm\ for r > a. Each mode is labeled as two mode num¬ 
bers (N, m) which are the radial and azimuthal mode 
numbers. The mode numbers N and m are also related 
to the number of maxima in the field profile in the ra¬ 
dial and azimuthal direction, respectively. Therefore, 
the field profile in Fig. [ 2 ] (c) corresponds to the mode 
(1, 2). From field intensity and Poynting vector distri¬ 
bution, we can conclude that the excitation of resonant 
particle modes lead to a negligible SPP transmission. In 
this case, the coupling of incident SPP to the WGMs 
of the particle is very similar to the critical coupling in 
the micro cavity-waveguide systems, in which the directly 
transmitted field and the field coupled out of resonator 
into the output waveguide are out of phase and thus van¬ 
ish each other [35] . We remark that, similar plasmonic 
transmission behavior has been reported in the SPPs 
scattering by deep narrow grooves structured into flat 
gold surfaces [21], for which efficient coupling of incident 
SPPs to resonant cavity modes inside the grooves leads 
to high reflectivity. It is also clear from Poynting vector 
pattern in Fig. [ 2 ] (c) that most incident power scatters 
to the radiation fields into the forward direction. 

In a simple geometric optics picture, modes of order 
m correspond to a plane wave with propagation constant 
of k = nko , undergo multiple total internal reflection 
upon the interior surface of the resonator. After one full 
circulation within the cavity (fi = 2tt in (J8|) we can write: 

kL rn + unfit = 27r m, (9) 

where L m = 2masin(^), is the feedback path and <fi t is 
the polarization-dependent phase shift that occurs dur¬ 
ing each total internal reflection [36]. Note that, the ray 
optics approximation is more suitable for large cavities, 
i. e. a A, however it can also describe low-order 
WGMs [22[. For a fixed wavelength we see from © that 
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L m oc m which show that higher mode number needs 
longer paths and therefore larger cavities. In fact, it is 
possible to excite higher order WGMs with increasing 
the particle diameter as demonstrated in Figs. [ 2 ] (e) and 
(g) (for (N,m)=(l, 3) and (1, 4)). The Poynting vector 
field also demonstrate how the field scattered in different 
directions for different particles. Figures [ 2 ] (d) and (f) il¬ 
lustrate the field patterns for nonresonant configurations 
where as can be seen most input power conveys to the 
transmitted surface plasmon polariton. So, the coupling 
of SPP mode to inter cavity modes of the particle leads 
dramatically changes in the transmission, reflection and 
scattering coefficients. This sensitivity of the transmit¬ 
ted SPPs to the particle size might be useful for single 
particle sensing down to subwavelength scales. In fact, 
the penetration of the evanescent near-field to the up¬ 
per dielectric medium makes the SPPs sensitive to the 
refractive index changes in their local environment. The 
presence of the particle can be easily determined because 
coupling between the SPP and particle causes a strong 
transmission change. Here, all-optical excitation and de¬ 
tection is one of the most important advantages of SPP 
sensing. 

We now investigate the optical force on the particles 
near a metal surface. Large field enhancement near 
a metallic surface allows efficient optical trapping us¬ 
ing evanescent plasmonic fields. Conventional optical 
traps have been widely used for manipulation of small 
dielectric particles and living cells, although the size of 
trapped objects are usually larger than the light wave¬ 
length. Plasmonic optical traps make it possible to ef¬ 
fectively trap smaller subwavelength particles on a flat 
surface at low input power [18]. Also, the surface plas¬ 
mon radiation pressure can be used for particle trans¬ 
port in a liquid environment Here, we show how the 
particle sizes strongly affect the scattering and gradient 
plasmonic forces. 

The normal and lateral optical forces has been demon¬ 
strated in Fig. [3] Again, we have two separate regains: 
small particle region where both F x and F z magnitudes 
are increased with increasing particle size and fluctuat¬ 
ing force region. However, for very small particles, due 
to strong field gradient of the SPP wave perpendicular to 
the surface the attractive normal force, F z is much larger 
than scattered lateral force, F x . In the fluctuating region 
we see the excitation of WG modes strongly affects the 
force magnitudes. We should note that similar behavior 
was reported from an incident plane wave on spherical 
dielectric particles [38]. However, for SPP incident wave 
the force fluctuations are much stronger compared to the 
plane wave incident. 

Next, we shortly explore the influence of particle shape 
to the transmittance and optical forces. The results are 
shown in Fig. [4] For comparison we again show the 
circular particle transmittance and forces. The size of 
square and equilateral triangle particle are define as the 
diameter of the circle circumscribing the particle. As 
can be seen, the particle shape has considerable influ- 



FIG. 3. (a) Lateral optical force, F x , and (b) normal optical 
force, F z , on a circular silicon cylinder as a function of its 
diameter near a flat silver surface at A = 600 nm. The gap 
between the particle and silver surface is 12 nm. 


ence on the transmittance. However, notice that, for 
some particle sizes (see, e.g., the arrow in Fig. 0(a)), 
all three particles show a resonant behavior, but the 
excited modes are different. Their magnetic field pat¬ 
terns are illustrated in Fig. |(e) and Figs jl] (b)-(c). It 
is possible to excite a standing wave mode for square 
particles. In fact, careful monitoring of the instanta¬ 
neous magnetic field over one period reveals that the 
mode profile of the particle of Fig. [4] (c) is a standing 
wave mode. The standing wave pattern can be given 
as H y (x',y') = A sin(m x / ttx'/o) sin(ra y /7n//a), where x' 
and y' are coordinates in the directions of the square sides 
of length a. We can thus deduce that Fig. [4] (c) shows 
the excitation of the mode ( m x /, m y >) = (2, 2). 

The induced optical forces are plotted in Figs. [4] (d)-(e) 
which again show high shape dependence. Particularly, 
F z is much bigger for circular particles compared to the 
square and triangular particles. 


V. CONCLUSIONS 

In conclusion, surface plasmons scattering by subwave¬ 
length defects and particles has been investigated using 
two dimensional finite difference time domain method. It 
is shown that even for small subwavelength particles the 
excitation of WG modes or standing wave modes strongly 
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affects the transmittance and reflectance of SPPs, which 
might have potential sensing applications. It is also 
shown that this mode excitation strongly influences the 
optical forces acting on the particle. While in this work 
we consider the 2D scattering problem, it can be straight¬ 
forwardly extended to more general structures, e.g., SPPs 
scattering by nanorods at oblique incidence. 


FIG. 4. (a) Transmittance as a function of particle size, (b)- 
(c) magnetic field patterns, \H y \ for the particles are shown 
by arrow in Fig. 4 (a), (d) the optical forces F x and (e) F z 
for three different particle shapes near a silver-water surface 
(A = 600 nm ). The gap between the particle (silicon) and 
silver surface is 12 nm. 
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